Renormalized interactions with a realistic single particle basis 
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Neutron-rich isotopes in the sdpf space with Z < 14 require modifications to derived effective 
interactions to agree with experimental data away from stability. A quantitative justification is given 
for these modifications due to the weakly bound nature of model space orbits via a procedure using 
realistic radial wavefunctions and realistic N N interactions. The long tail of the radial wavefunction 
for loosely bound single particle orbits causes a reduction in the size of matrix elements involving 
those orbits, most notably for pairing matrix elements, resulting in a more condensed level spacing 
in shell model calculations. Example calculations are shown for 36 Si and 38 Si. 
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I. INTRODUCTION 

New facilities for rare isotope beams will push the ex- 
perimental capabilities of nuclear physics with radioac- 
tive beams to more unstable, shorter-lived nuclei. Prop- 
erties of these nuclei exhibiting different behavior than 
stable nuclei, like the evolution of shell structure, are of 
significant interest for the next decades of research. A 
new theoretical technique and its behavior for stable and 
exotic nuclei has been studied to examine the importance 
of refining theoretical approaches for the production of 
model space interactions for unstable nuclei. 

Much research has been done using renormalization 
methods to convert a realistic interaction fit to nucleon- 
nucleon (NN) scattering data into an interaction in the 
nuclear medium. The goal is to renormalize the interac- 
tion to valence orbits outside of a stable, semi-magic or 
doubly magic nucleus treated as a vacuum in further cal- 
culations. A typical example would use 16 as the core 
and renormalize the NN interaction into the sd model 
space. For such an application, the harmonic oscillator 
basis of the form ^ n i mi (r) = -R^°(r) Fzm, (#,</>) is gener- 
ally used. Additionally, all the valence orbits are bound 
in the harmonic oscillator basis. For more exotic closed- 
subshell nuclei, loosely bound orbits often play a role. 
The harmonic oscillator basis is less applicable further 
from stability. Loosely bound orbits particularly deviate 
from the oscillator basis, as they exhibit a "long-tail" be- 
havior with a larger spread in the radial wavefunctions. 
However, few calculations have been done with a realis- 
tic radial basis for unstable nuclei with renormalized NN 
interactions. 

Experimental interest in neutron-rich silicon isotopes 
and the failure of some shell model Hamiltonians to re- 
produce data in the region have led to modifications in 
the SDPF-NR interaction pj, which had been the stan- 
dard for shell model calculations in the sdpf model space. 
The new SDPF-U interaction has different neutron- 
neutron pairing matrix elements for Z > 15 and Z < 14 



to account for the behavior of pf neutron orbits relative 
to the number of valence protons. The Z < 14 version of 
the interaction treats neutron-rich unstable nuclei that 
exhibit different shell behavior than the less exotic nu- 
clei in the Z > 15 nuclei. The interest in silicon isotopes 
and the nature of the SDPF-U interaction make 34 Si a 
suitable choice for the renormalization procedure with a 
realistic basis. A similar effect occurs for the neutron- 
rich carbon isotopes around the N = 14 closed subshell, 
requiring a 25% reduction in the neutron- neutron two- 
body matrix elements from the effective interactions de- 
rived for the oxygen isotopes Q ■ 



II. RENORMALIZATION PROCEDURE 

We begin with the realistic charge-dependent NN in- 
teraction N 3 LO derived at fourth order of chiral pertur- 
bation theory with a 500 MeV cutoff and fit to experi- 
mental NN scattering data [3J. The N 3 LO interaction 
is renormalized using a similarity transformation in mo- 
mentum space with a sharp cutoff of A = 2.2 fm _1 to 
obtain the relevant low momentum interaction We 
will refer to this technique as a vi ow k renormalization. 
Skyrme Hartree-Fock calculations are performed with the 
Skxtb interaction [5j for a chosen closed sub-shell target 
nucleus to determine the binding energy, single particle 
radial wavefunctions, and single particle energy spectra 
for neutrons and protons of the target nucleus. The low 
momentum interaction is then renormalized into a model 
space of interest using Rayleigh-Schrddinger perturba- 
tion theory [6] to second order including excitations up 
to 6huj, summing over folded diagrams to infinite order. 
We will compare three options for the renormalization 
to produce an effective interaction: harmonic oscillator 
single particle energies and wavefunctions (HO), Skyrme 
Hartree-Fock single particle energies and wavefunctions 
(SHF), and Skyrme Hartree-Fock single particle ener- 
gies and harmonic oscillator single particle wavefunctions 
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TABLE I: Single-particle energies for 34 Si and 40 Ca using the 
Skxtb interaction. Values in bold are in the model space. 
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(CP). 

The CP basis and HO basis give identical results to 
first order in perturbation theory since they use identi- 
cal wavefunctions. The energies, which are different in 
the two procedures, come into higher order diagrams via 
energy denominators, as discussed in (j . Therefore, the 
last option is a core-polarization basis (CP basis) since 
the core-polarization diagrams are affected to leading or- 
der even though the result at first order is unchanged. 

Skyrme Hartree-Fock radial wavefunctions, once 
solved, are implemented in the renormalization by using 
an expansion in terms of the harmonic oscillator basis 
via: 

^fV) = E^^MtW^ ® Xs]j, (1) 
n 

where a 2 gives the percentage of a specific harmonic os- 
cillator wavefunction component in the Skyrme Hartree- 
Fock wavefunction. The Skyrme Hartree-Fock wavefunc- 
tions and single particle energies can only be determined 
for bound states. For unbound orbits, the harmonic os- 
cillator basis remains in use, but the Gram-Schmidt pro- 
cess is used to ensure orthonormality of the single par- 
ticle wavefunctions. The effective interaction, consisting 
of the derived two-body matrix elements and the Skyrme 
Hartree-Fock single particle energies, can then be used in 
a shell model program directly. 

III. APPLICATION TO SDPF MODEL SPACE 

Neutron-rich silicon isotopes present an interesting ap- 
plication of the procedure outlined in the last section. A 
deeper understanding of the need for multiple interac- 
tions in the sdpf model space, as seen by the form of 
SDPF-U, can be gained by performing the renormaliza- 
tion for the same model space in multiple ways. The 



TABLE II: Single-particle energies for 34 Si and 40 Ca in the 
harmonic oscillator basis. The energy shift is chosen so that 
the valence energy is identical in both bases. Values in bold 
are in the model space. 
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model space chosen is the sd proton orbits and pf neu- 
tron orbits. The renormalization procedure is done using 
all three options for two different target nuclei, produc- 
ing a total of six interactions. The two target nuclei cho- 
sen are the stable 40 Ca doubly magic nucleus, and the 
neutron-rich 34 Si semi-magic nucleus. Single particle en- 
ergies of the SHF basis, using the Skxtb interaction, are 
presented in Table I for both target nuclei. For an SHF 
state that is unbound, the radial wavefunction is approx- 
imated by a state bound by 200 keV that is obtained by 
multiplying the SHF central potential by a factor larger 
than unity. The energy of the unbound state is estimated 
by taking the expectation value of this bound state wave- 
function in the original SHF potential. 

In the SHF basis, the calculation of single particle en- 
ergies shows that the proton orbits are shifted down in 
energy for 34 Si relative to 40 Ca, while the neutron or- 
bits are shifted up. For the valence neutrons, this shift 
results in a switch from four orbits for Ca bound by 
5.4 Mev on average to four orbits for 34 Si centered at 
0.0 MeV. This change, specifically the loosely bound en- 
ergies of the P3/2 and P1/2, has a significant effect on 
the wavefunctions, which will be discussed in more detail 
later. For comparison, the single particle energies used 
in the HO basis are given in Table II. The Blomqvist- 
Molinari formula hto = (45A- 1 / 3 - 25A~ 2 / 3 ) MeV 
gives 11.508 MeV for A =34 and 11.021 MeV for A = 
40. The absolute value of the harmonic oscillator basis is 
irrelevant, as only energy differences come into the dia- 
grams in Rayleigh-Schrodinger perturbation theory. For 
a better comparison to the SHF basis, the absolute value 
is chosen separately for protons and neutrons such that 

n va i 

(2 J+l)e a is identical in the HO and SHF bases, where 

1 

n va i, the number of valence orbits, is three for protons 
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40 Ca HO basis 

FIG. 1: Comparison of pf neutron- neutron matrix elements 
(in MeV) for the renormalization procedure in the HO basis 
for the two target nuclei. The solid line y = x denotes where 
the matrix elements would be identical. 
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40 Ca SHF basis 

FIG. 2: Comparison of pf neutron-neutron matrix elements 
(in MeV) for the renormalization procedure in the SHF ba- 
sis for the two target nuclei. The solid line y = x denotes 
where the matrix elements would be identical. Black dots 
correspond to matrix elements with J > 0, while the J = 
matrix elements are split into three groups: ff — ff (crosses), 
// — PP (diamonds), and pp — pp (plus signs). 



and four for neutrons and e a is the energy of the single 
particle orbit given by the a = n, l,j quantum numbers. 
In order to avoid divergences from the calculation of en- 
ergy denominators, all model space orbits are set to the 
same valence energy such that the starting energy Q of 
each diagram is constant. 



34 Si neutron radial wavefunctions 
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FIG. 3: Comparison of the single particle radial wavefunctions 
for 34 Si in the HO and SHF bases. 



Fig. 1 shows a comparison of the pf matrix elements 
in MeV for both target nuclei with the HO basis used 
in the renormalization procedure. The values deviate 
slightly from the line of equality but agree well with each 
other. Therefore, the choice of target nucleus, whether 
34 Si or 40 Ca, has little effect on the matrix elements in 
the HO basis. However, when we look at Fig. 2, where 
the comparison is for both target nuclei in the SHF basis, 
we see a reduction in the strength of the interaction for 
34 Si. This reduction with 34 Si as the target nucleus is due 
to the weakly bound nature of the pf neutron orbits. 

In the SHF basis, the / 7 / 2 orbit is bound by 2.6 MeV, 
and its radial wavefunction agrees well with the harmonic 
oscillator wavefunction as seen in Fig. 3. The Skyrme 
wavefunction is expanded in the harmonic oscillator ba- 
sis up to n — n max and the a n coefficients are renormal- 

ized to ensure that a„ = 1. For our renormalization 

71=0 

procedure, orbits up to (2n + I) = 9 arc included, which 
gives n max = 3 for the / 7/2 and / 5 / 2 orbits and n max = 4 
for the P3/2 and p\/2 orbits. This includes over 99% of 
the strength for the / orbits, but only 93% and 92% for 
the P3/2 and p\/2 orbits respectively. A first order calcu- 
lation can be done to n max = 6 for all orbits, which gives 
100%, 98%, and 97% for the /, P3/2, &ndpi/2 expansions 
respectively. 

With this procedure, 99% of the f 7 / 2 orbit is repre- 
sented by the R^f wavefunction. The pz/2 an d P1/2 
orbits are only bound by 400 and 269 keV, respectively. 
The expected harmonic oscillator component Ryf only 
makes up 80% and 78% of the respective radial wave- 
functions. Higher n orbits which extend farther away 
from the center of the nucleus contribute the remaining 
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TABLE III: Core polarization and total matrix elements in 
MeV of the form (aa\V\bb) j = o for different renormalization 
procedures. 
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strength. The / 5 / 2 orbit is unbound by three MeV, but 
the solution for the Skyrme radial wavefunction is deter- 
mined by assuming that the orbit is bound by 200 keV. 
With this method, 97% of the realistic radial wavefunc- 
tion is given by the R§f wavefunction. Single particle 
radial wavefunctions of valence space neutron orbits are 
shown in Fig. 3 in both the HO and SHF basis. The long 
tail behavior of the loosely bound p orbits can be seen 
in the SHF basis, as the wavefunctions have significant 
strength beyond 8 fm unlike the oscillator wavefunctions. 

The J = matrix elements in Fig. 2 deviate more 
from the line of equality, i.e. the pairing matrix ele- 
ments are reduced for 34 Si when the N 3 LO interaction 
is renormalized in the SHF basis. The SDPF-U inter- 
action has different neutron-neutron pairing matrix ele- 
ments for Z > 15 and for Z < 14 to account for 2p-2h ex- 
citations of the core correctly, depending on whether 34 Si 
or 40 Ca should be considered the core [J. The SDPF-U 
neutron-neutron pairing matrix elements are reduced by 
300 keV for Z < 14 in order to produce results in better 
agreement with experimental data. The pairing matrix 
elements in Fig. 2 are reduced for the 34 Si target by 213 
keV on average, relative to the case with 40 Ca as the tar- 
get. While the connection here to the Z-dependence in 
SDPF-U is only suggestive, the change in target mim- 
ics the change in core for calculations in the sdpf region, 
cited by Nowacki and Poves as the cause of their 300 keV 
reduction [![■ The reduction of 213 keV is due solely to 
the change in occupation of the d 5 / 2 proton orbit, which 
can affect the single particle energies and radial wave- 
functions, as well as the available diagrams in the core 
polarization. We find that the change in single particle 
radial wavefunctions plays the most significant role, but 
are also able to analyze the effect of the core polarization. 

Table III isolates a few matrix elements and compares 
the total matrix elements and the component due to core 
polarization for both target nuclei in all three bases. The 
reduction for total matrix elements involving p orbits is 
dramatic (« 30%) and is primarily due to the extension 
of wavefunction strength to large distances. Kuo et al. 
Q noted a reduction of core polarization in the harmonic 
oscillator basis and used different oscillator parameters to 



account for the core nucleons and valence nucleons sep- 
arately in halo nuclei. While 34 Si is not a halo nucleus, 
the loosely bound p orbits behave in much the same way 
as the valence nucleons in a halo nucleus. The reduction 
in core polarization is seen going from the 40 Ca target 
to the 34 Si target in any basis in Table III, although the 
size of the polarization is reduced for nucleons far from 
the core. As noted in the core interacts less with 
nucleons far away, so the excitations of the core are re- 
duced. The core polarization for matrix elements solely 
involving p orbits is under 100 keV. We observe that the 
core polarization can be reduced significantly without the 
total matrix element changing in the same proportion. 
For instance, the (f 7/2/7/2 [^1/7/2/7/2) matrix element is 
only reduced by 5% from 40 Ca to 34 Si in the HO ba- 
sis even though the core polarization is reduced by 30%. 
In the SHF basis, which takes into account the realis- 
tic wavefunction, the total matrix element is reduced by 
13% even though the core polarization is reduced by 43%. 
We would prefer to compare matrix elements involving 
the P3/2 or P1/2 orbits, but the core polarization becomes 
very small for loosely bound orbits, skewing percentage 
comparisons. Ogawa et al. [9] produce results which 
seem to be consistent with ours, identifying a 10%-30% 
reduction in nuclear interaction matrix elements involv- 
ing loosely bound orbits using a realistic Woods-Saxon 
basis. However, they were limited to comparisons of ra- 
tios of matrix elements and did not include core polar- 
ization. We show that core polarization suppression and 
reduction due to spread of the wavefunctions are both 
important effects which should be included, but do not 
tell the entire story. The / 7 / 2 wavefunction is very similar 
in the HO and SHF bases, as seen in Fig. 3, and yet the 
(/7/2/V/2 ^1/7/2/7/2) matrix element does not follow the 
same trend as the core polarization component in Table 
III. Other diagrams which are included at second order 
are relevant, and the full treatment of the renormaliza- 
tion in a realistic basis, as developed here, is necessary for 
accurate results. Our improvements enable us to perform 
calculations for neutron-rich silicon isotopes directly. 



IV. CALCULATIONS FOR 36 SI AND 38 SI 

The effect of the different interactions on nuclear struc- 
ture calculations has been studied as neutrons are added 
to Si. In order to obtain a consistent starting point, 
the proton-proton and proton-neutron matrix elements 
of SDPF-U have been used, with proton single particle 
energies (SPEs) chosen to reproduce those obtained by 
SDPF-U. Because SDPF-U does not reproduce the bind- 
ing energy of 35 Si, the SDPF-U neutron SPEs have been 
increased by 660 keV. The six interactions use neutron 
SPEs that reproduce the values of this modified SDPF-U 
interaction. 

The only difference in the six interactions used in the 
calculations are the neutron- neutron matrix elements. 
Calculations have been done in the model space discussed 
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FIG. 4: Calculations for the lowest energy states for J = 
0, 2, 4, 6 in 36 Si relative to 34 Si from the renormalization pro- 
cedure for 34 Si and 40 Ca, in the HO, CP, and SHF bases for 
both target nuclei. Crosses are used for calculations with 40 Ca 
as the target nucleus. 



in the last section with the shell model code NuShellx 
[i"oj . Fig. 4 shows the lowest J = 0,2,4,6 states in 
36 Si, relative to 34 Si. The HO basis and the CP basis 
for the same target nucleus deviate by no more than 20 
keV. However, the SHF basis noticeably shifts the states, 
with the largest effect being 170 keV more binding in the 
ground state with 40 Ca as the target nucleus. The bind- 
ing energy of 36 Si changes by nearly 500 keV depending 
on which renormalization procedure is used. Further- 
more, the level schemes for 36 Si are more spread out for 
the crosses where 40 Ca is chosen as the target nucleus. 

Fig. 5 shows the same states in 36 Si relative to 34 Si, 
but now the comparison includes the SDPF-U calcula- 
tions and experimental data. The CP basis results are 
not included since they are so similar to the HO basis cal- 
culations. We see that the level scheme for 36 Si is more 
spread out for the Z > 15 SDPF-U calculation than for 
the Z < 14 calculation, in agreement with our results 
discussed above. Our calculations for each method are 
in reasonable agreement with the comparable SDPF-U 
calculation, except for the + state which differs by over 
300 keV in both instances. The experimental binding 
energy relative to 3 Si is taken from a new mass mea- 
surement of 36 Si which is 140 keV higher in energy than 
previously measured (Tl| . The excitation energies of the 
Z < 14 SDPF-U calculation are comparable to experi- 



FIG. 5: Calculations for the lowest energy states for J = 
0, 2, 4, 6 in 36 Si relative to 34 Si using neutron-neutron matrix 
elements from SDPF-U and the renormalization procedure for 
both 34 Si and 40 Ca as target nuclei, using the SHF and HO 
bases. Experimental data is shown for comparison, with a 
new mass from [TJ]. Crosses are used for calculations with 
40 Ca as the target nucleus. 



ment, as expected from an interaction fit specifically to 
neutron-rich silicon isotopes. While no interaction repro- 
duces the experimental data very well, general trends can 
be seen. The calculations with 40 Ca as the target nucleus 
depicted by crosses result in level schemes that are too 
spread out in comparison to the experimental data. The 
reduction in the strength of the interaction for 3 Si using 
the SHF basis results in a reduction of the energy of the 
states in 36 Si, especially for the ground state (the pairing 
matrix elements were most reduced). The rms between 
experiment and theory with 34 Si as the target nucleus in 
the SHF basis is 214 keV for the four states shown. One 
reason for this deviation is the lack of three body forces 
in the procedure. The inclusion of the NNN interaction, 
at least via the effective two-body part, is important for 
a higher level of accuracy. Additionally, the chosen SPEs 
may contribute to the deviation, which would be better 
constrained if all the single particle states in 35 Si were 
known experimentally. For exotic nuclei, the calculated 
single particle state is often above the neutron separa- 
tion energy and determination of the experimental single 
particle states may not be possible with current facilities. 
Thus it is essential to improve energy density functionals 
such that they provide reliable single particle energies. 
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FIG. 6: Calculations for the lowest energy states for J = and 
J = 2 in 38 Si relative to 34 Si using neutron-neutron matrix 
elements from SDPF-U and the renormalization procedure 
for both 34 Si and 40 Ca as target nuclei, using the SHF and 
HO bases. Crosses are used for calculations with 40 Ca as the 
target nucleus. 

The importance of using a realistic basis and an appro- 
priate target nucleus to renormalize an interaction into 
the nuclear medium for calculations of exotic nuclei is ev- 
ident from the calculations shown here. Otherwise, the 
interaction will be too strong and will produce signifi- 
cant overbinding even for the two particle case, an effect 



that gets magnified as more particles are added, as seen 
in Fig. 6 for 38 Si. Only the + and 2 + states are shown 
since the 4 + and 6 + states are not known experimentally, 
but the binding energy is only approximately reproduced 
for the calculations with 34 Si in the SHF basis. As noted 
in the 36 Si case, the excitation energy of the 2 + state is 
too high in the SHF basis but is reproduced well by the 
Z < 14 SDPF-U calculation for 38 Si. 



V. SUMMARY AND CONCLUSIONS 

The microscopic nucleon-nucleon interaction N 3 LO 
was renormalized using vi ow k and many-body perturba- 
tion theory in order to produce an effective interaction in 
the nuclear medium that could be used in a shell model 
code. The renormalization was performed in three dif- 
ferent bases: harmonic oscillator, core polarization, and 
Skyrme Hartree-Fock. The choice of basis can signifi- 
cantly affect the value of matrix elements, as shown in 
the comparisons of pf neutron-neutron matrix elements 
for the stable 40 Ca and the neutron-rich 34 Si nuclei. The 
difference results from valence orbits being bound by only 
a few hundred keV, resulting in a long tail in the radial 
wavefunction relative to the harmonic oscillator wave- 
function. The loosely bound orbits cause a reduction 
in the overall strength of the interaction, an effect that 
becomes magnified as full scale shell model calculations 
are performed. Accounting for the properties of the or- 
bits by using a realistic basis is essential for an accurate 
description of the nuclear interaction in exotic nuclei as 
determined by the renormalization of an NN interaction. 
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